Andreev Current-Induced Dissipation in a Hybrid Superconducting Tunnel Junction 
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We have studied hybrid superconducting micro-coolers made of a double Superconductor- 
Insulator-Normal metal tunnel junction. Under subgap conditions, the Andreev current is found 
to dominate the single-particle tunnel current. We show that the Andreev current introduces addi- 
tional dissipation in the normal metal equivalent to Joule heating. By analyzing quantitatively the 
heat balance in the system, we provide a full description of the evolution of the electronic temper- 
ature with the voltage. The dissipation induced by the Andreev current is found to dominate the 
quasiparticle tunneling-based cooling over a large bias range. 

PACS numbers: 74.50.+r, 74.45.+C 
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In a tunnel junction between a normal metal (N) and 
a superconductor (S), the charge transfer occurs mainly 
through two different mechanisms. The tunneling of a 
single quasi-particle is possible for electrons or holes with 
an energy E (compared to the Fermi level Ep) larger than 
the superconductor gap A. At low energy, the charge 
transfer occurs through the Andreev reflection 

0,1- In 

the normal metal, an electron (a hole) impinging on the 
superconducting interface is reflected as a hole (an elec- 
tron), enabling the transfer of a Cooper pair into (out 
of) the superconductor. As the energies of the involved 
electron and hole are located symmetrically around Ep, 
the Andreev reflection does not carry heat through the 
interface at zero bias. The probability for an incident 
quasi-particle to follow an Andreev reflection, a specu- 
lar reflection or a tunnel transfer is given in the ballistic 
regime (no disorder) by the BTK theory [3J. For a N-I-S 
tunnel junction with an insulator (I) of intermediate or 
low transparency, the Andreev reflection probability is 
predicted to be vanishingly small. Taking into account 
the quasi-particles confinement in the vicinity of the in- 
terface, this is no longer true. This confinement can be 
induced by the disorder or the presence of a second bar- 
rier in the normal metal. A single quasiparticle then ex- 
periences several collisions with the interface 0, [5] • The 
actual Andreev reflection transmission coefficient corre- 
sponds to the coherent addition of many individual trans- 
mission probabilities. Therefore, the Andreev sub-gap 
current significantly exceeds the ballistic case prediction 
[H, 0] and can be modulated by a magnetic flux [j| . 

A quasi-particle current in a N-I-S junction indeed car- 
ries both a charge current and a heat current. With a 
voltage bias smaller than the gap A/e, the tunnel cur- 
rent is selectively made out of high-energy electrons (or 
holes) ; this cools the electronic population of the normal 
metal . In this way, (S-I- N-I-S) micro-coolers based on 
a double tunnel junction provide a significant tempera- 
ture reduction which reaches an optimum at a voltage 



bias just below the gap. At a very low temperature, the 
thermal transport in such N-I-S tunnel junctions appears 
to be still little understood. For instance, an apparent 
reversal of the normal metal temperature evolution was 
observed in various experiments 1(J ll[ and related to a 
non-BCS density of states of the superconductor A 
clear understanding of this behavior is still missing. 

In this letter, we describe an experimental study of the 
heat transport in a S-I-N-I-S junction, focusing on the 
very low temperature regime. We show that the phase- 
coherent Andreev current introduces a significant dissi- 
pation in the normal metal. We provide a fully quanti- 
tative analysis of the heat transfer in the system which 
shows that although the Andreev current is a small ef- 
fect in terms of charge current, the heat it creates has a 
dominating influence on the heat balance. 

Fig. Q] inset shows a typical sample which features 
a geometry similar to the one studied in Ref. [ijj]- It 
consists of a 50 nm thick, 4 /jm long and 0.3 /im wide 
normal metal Cu electrode embedded between two 40 nm 
thick and 1.5 /im wide superconducting Al electrodes. 
The tunnel barriers at the two symmetric junctions of 
dimensions 1.5 x 0.3 fim 2 were prepared by oxidization in 
0.2 mbar of O2 pressure for 3 min and give a total normal- 
state resistance R n in the range 2-3 kO. In addition to the 
two cooling junctions, we added three Cu tunnel probes 
of area 0.33x0.43 /im 2 on one Al electrode (one is shown 
in Fig. [1] insets). These probes are strongly connected 
to both a Cu reservoir and an Al reservoir so that no 
cooling is expected there. These junctions thus provide a 
reference for an isotherm behavior. In the following, we 
will describe the behavior of one out of three investigated 
samples which all showed similar behavior. 

We have measured the current-voltage I(V) charac- 
teristic of every probe junction and of the two cooling 
junctions in series at temperatures down to 90 mK. The 
differential conductance dI/dV(V) is obtained by numer- 
ical differentiation. We have taken special care to obtain 
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FIG. 1: (Color online) Left inset: Geometry of the sample. 
Right inset: Scanning electron microscope micrograph of a 
typical cooler sample made of a normal metal Cu electrode 
(light grey) connected to two superconducting Al reservoirs 
(dark grey) through tunnel junctions. One of the additional 
probe junctions connected to one Al reservoir is visible at the 
bottom. Main figure: normalized differential conductance as 
a function of the voltage and at the cryostat temperatures of 
90 (red curve), 230 (blue), 330 (green) and 440 mK (brown). 
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FIG. 2: (Color online) Current-voltage characteristic of the 
cooler junctions (blue symbols) and of the probe (red symbols) 
as a function of the voltage and at a cryostat temperature of 
90 mK together with best-fit calculations of the Andreev cur- 
rent (dotted lines) and of the single-particle current (dashed 
lines) using the parameters D = 80 cm 2 .s _1 , L v = 1.5 fim, 
R n = 1.9 kQ. (cooler) and 2.76 kQ (probe), 2A = 0.43 meV 
(cooler) and A = 0.228 meV (probe). The fit temperature is 
105 mK. Compared to the theoretical calculations based on 
Ref. the Andreev current was multiplied by 1.55 (cooler) 
and 4.3 (probe). No cooling effect is included in the calcula- 
tion. 



accurately the sub-gap conductance of our current-biased 
samples down to a level of about 10 -4 of the normal state 
conductance. Fig. Q] displays on a logarithmic scale the 
differential conductance of the cooling junctions for a se- 
ries of cryostat temperatures. The tunnel current in a 
N-I-S junction is given by: 

eR n Jo 2 2 

(1) 

where V/2 is the voltage across the N-I-S junction, Jn(E) 
is the electron energy distribution function in the nor- 
mal metal and n s (E) =| E | /%/ E 2 — A 2 is the normal- 
ized BCS density of states in the superconductor. An 
isotherm N-I-S junction would then feature in the sub- 
gap regime a linear behavior of the differential conduc- 
tance on a logarithmic scale. In contrast, we observe in 
the high temperature (T > 200 mK) data an upward cur- 
vature which constitutes a clear signature of the electron 
cooling [lij ]. 

In the low temperature regime (T < 200 mK), clearly 
a different characteristic is obtained with a differential 
conductance peak at zero bias. A similar behavior is 
obtained on every probe (not shown). This zero-bias 
peak cannot be accounted for by a single-particle tunnel 
current. The zero-bias differential conductance increases 
while the temperature is lowered below about 200 mK, 
which suggests that it is a phase-coherent effect. 

We will ascribe the low bias differential conductance 



peak to an Andreev current, i.e. a double particle tun- 
nel current created by Andreev reflections at the N-I-S 
junctions. In order to calculate the Andreev current I a, 
we used the theory of Ref. We took into account 

the finite gap A and the disorder both in the normal 
metal and in the superconductor. We considered the ID 
regime where the coherence length of an Andreev pair in 
the normal metal Le = \/TiD/E [13[ is much larger than 
the junction dimension. 

Let us first discuss the current-voltage characteristics 
of one probe junction, where the electronic tempera- 
ture can be considered as constant and very close to 
the cryostat temperature. Fig. [2] displays the mea- 
sured current-voltage characteristic of one probe junction 
(1.55 /im away from the cooler junctions) together with 
the calculated single quasi-particle and Andreev currents. 
We have fitted the current-voltage characteristics of the 
probe while taking the electron diffusion coefficient in 
Cu equal to the measured value D — 80 cm 2 .s _1 . The fit 
parameters are A = 0.228 meV, the electronic temper- 
ature T e — 105 mK, which is found slightly higher than 
the bath temperature, and L v = 1.5 fj/m, which agrees 
well with expectations for a pure metal at very low tem- 
perature [LJ]. In this fit, we had to scale the Andreev 
current by a multiplying factor M— 4.3 in a similar way 
to Ref. 8]. The fit describes very well the probe data. 
The cross-over between the low-bias Andreev current and 
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FIG. 3: (Color online) Current-voltage characteristic of the 
S-I-N-I-S cooler junction as a function the voltage at cryo- 
stat temperature of 90 (red line), 230 (blue) and 330 mK 
(green) together with best-fit calculated curves. Dot-dashed 
line: model including the cooling by the tunnel current, with 
parameters 2A = 0.43 meV, K.A = 144 W.K~ 4 . Dashed 
line: model including in addition the Andreev current, but 
not the related heating, with the parameters D — 80 cm 2 .s _1 , 
Lip =1.5 fim,M = 1.55. Dotted lines: full model taking into 
account the Andreev current and the related heating, 



the high-bias single quasi-particle current is clearly vis- 
ible. We did not need to take into account the possible 
contribution of pinholes in the barrier, which means that 
our junctions can be considered as ideal. We conclude 
that the sub-gap current in our N-I-S junctions is the 
superposition of a single-particle tunnel current and a 
two-particle Andreev current. 

We tried to fit the S-I-N-I-S cooler junction charac- 
teristic in a similar way. Fig. [2] also shows the current- 
voltage characteristic of the cooler junction together with 
the calculated single quasi-particle and Andreev currents. 
A good fit to the theory is found only at small bias, where 
cooling or heating effect remain negligible. This fit is ob- 
tained with the same parameter set than for the probe, 
except for the gap 2A = 0.43 meV and the scaling factor 
M = 1.55. The difference in the M factor is presum- 
ably due to difference in geometry between the types of 
junctions. 

In order to analyze quantitatively the behavior of our 
S-I-N-I-S cooler junctions, we need to consider the heat 
balance in the normal metal. Here we assume a quasi- 
cquilibrium situation: the electrons and the phonons fol- 
low a thermal energy distribution functions at a respec- 
tive temperature T e and T p h which are in general different 
from the cryostat temperature Tf, ase - This is justified be- 
cause the inelastic scattering time is of the order of the 
phase-coherence time of 300 ps (derived from Fig. 2 fit), 
which is much shorter than the mean escape time from 



the island estimated to about 100 ns. The tunnel current 
is responsible for a cooling power in the normal metal: 

Pcooi = 4w r (E-^)n s (E)[f N (E-^-)-f s (E)]dE 

e Rn J-oc 2 2 

(2) 

where fs is the energy distribution function in the su- 
perconductor. This cooling power is compensated by the 
electron-phonon coupling power so that: 2P CO oi+Pe- P h = 
0. Here the factor 2 accounts for the presence of two 
junctions. We used the usual expression for the electron- 
phonon coupling P e - ph (T e , T ph ) = EJ7(T e 5 - T* h ), where 
S is a material-dependent constant and U is the metal 
volume. If one considers now the normal metal phonons, 
the electron-phonon coupling power is compensated by 
the Kapitza power P K (T ph ,T base ) = KA(T* ase - T p 4 J, 
where K is an interface-dependant parameter and A the 
contact area, so that P e - P h + Pk = 0. The related 
Kapitza thermal resistance is significant only at inter- 
mediate temperature and above (T > 300 mK), in which 
case the normal metal phonons can be cooled below the 
substrate temperature 12j |. 

Fig. [3] displays a series of current-voltage character- 
istics obtained at different temperatures together with 
calculated curves. We have first numerically solved the 
heat balance equations in the high temperature regime 
where the Andreev current contribution can be neglected. 
We took the well-accepted value £ = 2 nW./im _3 .K -5 
[l5| and we determined from the fitting procedure the 
Kapitza coupling parameter value K.A — 144 W.K~ 4 . 

Fig. [3]dash-dotted line displays the calculated current- 
voltage characteristic calculated at a 90 mK cryostat 
temperature by including the charge and heat currents of 
only the single-particle channel effects. The agreement is 
poor, which confirms the need for including the Andreev 
current channel. As a second step, we included the An- 
dreev charge current component. The dashed line shows 
the result of the calculation with again the same param- 
eters set than in Fig. [2] and a cryostat temperature of 90 
mK. It provides an acceptable fit at low bias but shows 
a clear discrepancy at intermediate voltage. Adding a 
leakage term in the heat balance with a linear resistance 
does not provide a good description of the data. Thus 
a significant thermal contribution is missing in the heat 
balance equations described above. 

We investigated theoretic ally the heat transfer created 
by the Andreev current I a [16( flowing through an N-I-S 
tunnel junction (l7| . straightforwardly extending We 
have found that the work performed by the current source 
feeding the circuit with a current I a generates a Joule 
heating IaV that is deposited in the normal metal. The 
full heat balance equation for the normal metal electrons 
can then be written as: 

ZPcooi + Pe- P h + IaV = 0. (3) 
With this complete heat balance equation taken into ac- 
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FIG. 4: (Color online) Dependence of the calculated electronic 
temperature with the voltage as derived from the fit of the 
experimental result and for a series of cryostat temperature: 
90 (red symbols), 140 (green) and 230 mK (blue). 



count, we calculated the current- voltage characteristic at 
every cryostat temperature, see the dotted lines fits in 
Fig. [3] The agreement is very good for every accessible 
voltage and temperature. Assuming a thermalization of 
the phonons to the substrate temperature would change 
the current by less than 2% at 90 mK, which means that 
the phonon cooling has a negligible role in the data analy- 
sis presented here. Let us note that at the 90 mK cryostat 
temperature, the agreement covers 4 orders of magnitude 
for the current. 

This excellent agreement demonstrates that the heat 
balance in our S-FN-I-S junctions can be fully under- 
stood by taking into account the contributions of both 
the single quasi-particle tunneling and the two-particle 
Andreev current. The numerical solution of the heat 
balance equations also provides the electron temperature 
at every bias. Fig. 0] shows the calculated dependence 
of the normal metal central island electron temperature 
with the voltage across the cooler. At very low temper- 
ature, the electron temperature first increases with the 
bias due to Andreev current heating, before decreasing 
due to the tunnel current-based cooling effect. At a cryo- 
stat temperature of 90 mK, the cooling effect overcomes 
the Andreev current thermal contribution only close to 
the optimum bias. This demonstrates the great impor- 
tance of the Andreev current-induced dissipation in a N- 
I-S tunnel junction. Although the Andreev current is a 
small effect in a such junction if one considers the charge 
current, this is no longer true if one considers the heat 
current. The explanation is the following. Compared to 
a Joule power I.V, the Andreev current contributes fully 
to heating while the tunnel current cools with a moder- 
ate efficiency This efficiency is of the order of T e /A 
which is about 5 % at a 100 mK electron temperature. 



In conclusion, we have devised a quantitative analysis 
of the current- voltage characteristics of S-I-N-I-S hybrid 
junctions at very low temperature. Our study demon- 
strates the importance in terms of heat transport of the 
Andreev current which arises from the confinement of 
phase-coherent quasi-particles in the vicinity of the inter- 
face. This Andreev current-induced dissipation is with no 
doubt also of great importance in the behavior of S-N-S 
junctions with transparent interfaces, where a significant 
electron heating can be observed [l8| . 
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